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Abstract 

We prove collinear factorization theorem for the process irj* — > tt at the twist-3 level in the co- 



■ variant gauge by means of the Ward identity, concentrating on the two-parton case. It is shown that 

>. 

soft divergences cancel and collinear divergences are grouped into the pseudo-scalar and pseudo- 

(N 

tensor two-parton twist-3 pion distribution amplitudes. The delicate summation of a complete 

<n : 

© . S et of diagrams for achieving factorization in momentum, spin, and color spaces is emphasized. 

The proof is then extended to the exclusive semileptonic decay B — ► nlu, assuming the hard scale 
"OhI to be of 0(VAM B ), where A is a hadronic scale and Mb the B meson mass. We explain the 

distinction between the factorization of collinear divergences for a pion distribution amplitude and 



of soft divergences for a B meson distribution amplitude. The gauge invariance and universality 



of the two-parton twist-3 pion distribution amplitudes are confirmed. The proof presented here 
can accommodate the leading twist-2 case. We then compare our proof with that performed in the 
framework of soft-collinear effective theory. 

PACS numbers: 12.38.Bx 
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I. INTRODUCTION 



Recently, we have proposed a simple proof of collinear factorization theorem in pertur- 
bative QCD (PQCD) for the exclusive processes iry* — ► Tf7r) and B —> j(ir)W based on 
the Ward identity l|. According to this theorem S hadronic form factors are 

factorized into the convolution of hard amplitudes with hadron distribution amplitudes in 
momentum, spin, and color spaces. The former, being infrared finite, are calculable in per- 
turbation theory. The latter, absorbing the infrared divergences involved in the processes, 
are defined as matrix elements of nonlocal operators. The universality of the distribution 
amplitudes and the gauge invariance of the factorization have been explicitly demonstrated. 
Our proof can be compared to that performed in the axial gauge 3], in which the fac- 
torization of infrared divergences is trivial, but the gauge invariance is not obvious. The 
formalism in 1] is restricted to the leading- twist, i.e., twist-2 level. As emphasized in 0,0], 
contributions from the two-parton twist-3 pion distribution amplitudes are not only chi- 
rally enhanced, but of the same power as the leading-twist one in the semileptonic decay 
B — ► irlv. Hence, it is necessary to derive the corresponding factorization theorem. This 
proof can be regarded as an essential step toward a rigorous construction of factorization 
theorem for two-body nonleptonic B meson decays. 

The general decompositions of the matrix elements relevant to the two-parton pion dis- 
tribution amplitudes are, quoted from p, 

(0|%)7 M 7 5 n(0)|7r + (P)) = iUP, [ dxe~ ix p %(x) + kM^ f dxe'^g^x) (,1) 

Jo 1 " ■ y Jo 

(0|%) 7 5^0)|7r + (P)> = -if n m [ dxe- lxP - y (t>s{x) , (2) 

Jo 

i ( M 2 \ r 1 
(0|J( 2 /) 7 5^^(0)|7T + (P)) = --Umo 1 - HH (P*V» - PvVv) / dxe-™ p -y<j> a {x) , (3) 

" V m o / 

where <pv,s,a an d g n are the distribution amplitudes of unit normalization, f w the pion decay 
constant, M n the pion mass, x the momentum fraction associated with the d quark evaluated 
at the coordinate y. The Wilson links that render the above nonlocal matrix elements gauge 
invariant are not shown explicitly. It is easy to observe that the contribution from <pv, 
independent of the pion mass, is twist-2, and the contribution from g n is twist-4 because 
of the factor M%. The contributions from the pseudo-scalar (PS) distribution amplitude 
<ps and from the pseudo-tensor (PT) distribution amplitude <p a , proportional to the chiral 
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enhancing scale mo, are twist-3. 

We concentrate on the factorization of the two-parton twist-3 distribution amplitudes (fis 
and CT from the processes rcy* — > it and B — > ttW. We shall not consider the three-parton 
twist-3 distribution amplitudes here, since their contributions to exclusive processes are 
suppressed by the strong coupling constant and of higher power. The reason is as follows: 
after factorizing the corresponding infrared divergences to all orders, the contribution is 
written as a convolution of a hard amplitude with the three-parton twist-3 distribution 
amplitudes. The hard amplitude contains one more attachment from the extra parton 
(gluon) compared to that in the two-parton case. The attachment introduces one more 
power of the coupling constant, and one more hard propagator proportional to 1/Q, where 
Q is a large scale characterizing the hard amplitude. Hence, a three-parton hard amplitude 
is at least down by a power of the coupling constant and a power of 1/Q compared to the 
leading-order leading-twist hard amplitude. Moreover, the three-parton twist-3 distribution 
amplitudes should be considered along with the two-parton kx distribution amplitudes, 
which form a complete gauge-invariant set. 

Nonperturbative dynamics is reflected by infrared divergences of radiative corrections in 
perturbation theory. There are two types of infrared divergences, soft and collinear. Soft 
divergences come from the region of a loop momentum /, where all its components diminish. 
Collinear divergences are associated with a massless quark of momentum P ~ (Q, 0, 0^). In 
the soft region and in the collinear region with I parallel to P, the components of I behave 
like 

r = (/+, T, It) ~ (A, A, A) , l» ~ (Q, A 2 /Q, A) , (4) 

respectively, where the light-cone coordinates have been adopted, and A is a small hadronic 
scale. In both regions the invariant mass of the radiated gluon diminishes as A 2 , and the 
corresponding loop integrand may diverge as 1/A 4 . As the phase space for loop integration 
vanishes like dH ~ A 4 , logarithmic divergences are generated. 

In this paper we shall derive the collinear factorization formula for the scattering process 
7T7* — > 7r, which involves the pion form factor, at twist-3 by means of the Ward identity. 
The chirally enhanced contributions to the pion form factor have been calculated in [l_Q(] 
without proving their factorization theorem. It will be shown that soft divergences cancel 
and collinear divergences, factored out of the processes order by order, are absorbed into the 
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two-parton twist-3 pion distribution amplitudes denned by the nonlocal matrix elements, 



MX = - — 



P + 

1 d 
6 dx 
P 



dy 

2tt 



(0\d(y-) l5 V exp 



-w 



dzn^ ■ A(zri-) 



u(0)\tt + (P)) , 



x 



m 



dy 

2tt 



"~(0|d(y-) 75 (^ + tf_-l)Pexp 



dzn_ ■ A(zn-) 



where n + = (1, 0, 0-r) and n_ = (0, 1, 0t) are dimensionless vectors on the light cone, the 
symbol V stands for path-ordering of the Wilson line, and the pion decay constant f n has 
been omitted. The definition of the hard amplitudes at each order will be given as a result 
of the proof. 

We then prove the collinear factorization theorem for the semileptonic decay B — > ttW, 
whose topology is similar to the scattering process 7c-f* — > n. In the heavy quark limit the 
mass difference between the B meson and the b quark, A = Mb — rrib, represents a small 



scale. Assuming the hard scale to be of 0(a/AM#), the soft divergences do not cancel on 
the B meson side, and the B meson distribution amplitudes are introduced to absorb the 
soft divergences. The distinction between the factorization of soft divergences for the B 
meson distribution amplitudes and the factorization of collinear divergences for the pion 
distribution amplitudes will be explained. It will be shown that the two-parton twist-3 pion 
distribution amplitudes derived from the scattering 7T7* — > tt and from the decay B — > ttW 
are identical as defined by Eq. (j3J). That is, the universality of hadron distribution amplitudes 
is confirmed. 

There are different opinions on whether the transverse degrees of freedom of partons 
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The conclusion drawn in 



should be involved in exclusive B meson decays 

[12I Q] that the parton transverse momenta kr are not necessary is based on the analysis 
of the B — > decay, for which the collinear factorization formula does not develop an 

la ], for example, in the 
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end-point singularity. When end-point singularities appear [lj, 
collinear factorization formulas of semileptonic and nonleptonic decays, the region with a 
small momentum fraction x becomes important. In this end-point region the parton kx 
should be taken into account, and kr factorization theorem 12, Q] is more appropriate. 
Here we shall derive the collinear factorization formalism for exclusive B meson decays. The 
kx dependence can be introduced straightforwardly following the procedure in [lflj ]. 
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We emphasize that the proof of factorization theorem is not the whole story of PQCD. 



The double logarithms a s In x 
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appeari ng i n higher-order corrections to exclusive B meson 
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221 ] . When the end-point region is important, 



decays have been observed 
a s In 2 x can not be treated as a small expansion parameter, and should be summed to all 
orders. A systematic treatment of these logarithms has been proposed by grouping them 
into a quark jet function 



231 ] , whose dependence on x is governed by an evolution equation 
2l| . A Sudakov factor, obtained by solving the evolution equation, decreases fast at the end 
point. Moreover, if kj- factorization theorem is adopted, resummation is also required, 
which leads to another Sudakov factor describing the parton distribution in kr- Therefore, 
in a self-consistent analysis the original factorization formulas should be convoluted with 



;he above two Sudakov factors. It turns out that the end-point singu 



24 1 , and an arbitrary infrared cutoff for the momentum fraction x 



14. 



arities do not exist 



251 ] is not required. 



In the framework of soft-collinear effective theory (SCET) 2y,|27(, an effective Lagrangian 
with high-energy modes integrated out has been constructed. This SCET Lagrangian pro- 
vides a simple guideline for deriving a factorization formula by counting the powers of 
effective operators: start with an effective operator relevant for a high-energy QCD process, 
and draw the diagrams based on the SCET Lagrangians. Those effective diagrams, whose 
contributions scale like the power the same as of the operator, contribute to the matrix 
element formed by the operator. This matrix element is identified as the nonperturbative 
distribution amplitude, which collects the infrared divergences in the process. The Wilson 
coefficient (the hard amplitude) associated with the considered operator is then obtained 
by subtracting the effective diagrams from the full diagrams (the matching between the full 
theory and the effective theory). An example for the application of SCET, the collinear 
factorization of the B — > Dn decays, can be found in |28|. We shall compare the above 
formalism with ours at the end of this paper. For a detailed comparison, refer to j^. 

We mention the opinion from the QCD-improved factorization (QCDF) j^j], which claims 
that the B — > tt form factor, suffering the end-point singularity in collinear factorization, 
is dominated by soft dynamics. For a debate on this issue, refer to 0]. We have 
explained that the opposite conclusions on the dominant dynamics in exclusive B meson 
decays are attributed to the different theoretical frameworks j^lj]: a transition form factor 
is factorizable in PQCD, i.e., in k^ factorization as explained above, not factorizable in 
QCDF, i.e., in collinear factorization (speaking of only the leading contribution) j^j], and 
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partially factorizable in SCET |32| (non-singular and singular pieces in collinear factorization 
are written into factorizable and nonfactorizable forms, respectively). Hence, there is no 
conflict at all among the above observations. It has been pointed out that the collinear 
factorization and the kx factorization lead to different phenomenological predictions for 
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33, 



nonleptonic B meson decays, such as the CP asymmetries in the B — > tt + tt~ modes 
0]. Therefore, it is expected that the comparison with experimental data can discriminate 
the above approaches. 

In Sec. II we derive the 0(a s ) factorization of the collinear divergences in the process 
7T7* — > 7i. The delicate summation of a complete set of diagrams for achieving the factoriza- 
tions in momentum, spin, and color spaces is emphasized. The all-order proof based on the 
Ward identity is presented in Sec. III. The absence of the soft divergences is also shown. The 
technique is then generalized to the decay B — > ixlv in Sec. IV. The factorizations of soft 
divergences for the B meson distribution amplitudes and of collinear divergences for the pion 
distribution amplitudes are compared. Section V is the conclusion. We refer the detailed 
calculations of the 0(a s ) infrared divergences in the above two processes to Appendices A, 
B, and C. 



II. 0(a s ) FACTORIZATION OF vr 7 * -> vr 

We start with the two-parton twist-3 factorization of the process Tc-f* — > tt at the one-loop 
level, which will serve as a starting point of the all-order proof. The momentum Pi (P 2 ) of 
the initial-state (final-state) pion is parameterized as 

P x = (P+, 0, Or) = 0, Or) , P 2 = (0, P 2 , Or) = -%0, 1, T ) . (6) 

Consider the kinematic region with large Q 2 = —q 2 , q = P2 — P\ being the momentum 
transfer from the virtual photon, where PQCD is applicable. The lowest-order diagrams 
with the valence quarks being the external particles are displayed in Fig. 1. The lower 
valence quark (an anti-quark d) in the initial state pion carries the fractional momentum 
x\P\. The lower valence quark in the final state carries the fractional momentum £2-^2- 
Figure 1(a) gives the amplitude, 

n t$)( \_ i 2n di{xiP 1 )Ii j [Yd(x2P2)u{x 2 P2)'lv{p2 - Xi Pl)^]jlIlkUk(XlPl) 

Lr' '{x 1 ,x 2 — -eg U F — — — — — , {( 

2 (P 2 - XiPi) 2 (xiPi - X2P2Y 
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with xi( 2 ) = 1 — xi(2), the identity matrix I, and the color factor Cp, where the averages 
over spins and colors of the u and d quarks have been done. The u and d quark fields obey 
the equations of motion, 



xi fru{x x Pi) = , d(x 1 P 1 )x 1 A = , 

where the quark masses m u and have been neglected in high-energy scattering. 
Insert the Fierz identity, 



(8) 



Ii-ili 



^Iikhj + ^Cla)ik{i a )lj + ^(75 i-)ik{i+l$)lj 

+^(7s)ifc(75)ii + i- - i- ~ 1 )7s]^i 



(9) 



into Eq. (JJJ) to separate the fermion flow. Different terms in the above identity correspond 
to contributions of different twists. The PS structure proportional to 75 and the PT struc- 
ture proportional to 75(^1+ jh- — 1) lead to the twist-3 contributions. Equation (JJJJ) is a 
modified version appropriate for extracting collinear divergences (2^: the choice of the PT 
structure in Eq. © and the ordinary one ( , j5(J al3 )ik((Jai3l5)ij are equivalent. The PS and PT 
contributions to the process 717* — > 7r must be included simultaneously in order to form the 
gauge interaction vertex of a pseudo-scalar particle, which is proportional to (Pi + P^)^- 
The insertion on the initial-state side gives 

^\x 1 ,Ci)HP(C l ,x 2 )+^\x 1 ,Ci)H?(C 1 ,x 2 ) 



GV>(x lt 



X 2 



(10) 



The functions 



do) 

J S(T) 



and H 



(o) 

S(T)' 



<Pt\xi,^i) 



Ht\{,i,x 2 ) = l -eg 2 C F m 



1 



4mo 
1 



4m 



d{x 1 P l )^ h u{x x P 1 )8{$, 1 - xx) , 

J(xiPi)7 5 (^ + ^_ - 1)w(xiPi)5(£i - xi) , 

tr[ 7 ^(x 2 P 2 )M(x 2 P 2 )7,(P 2 -£1 A)7m7s] 



(P2 - 6 -Pi) 2 (6 -Pi - x 2 P 2 ) 2 
tr[Yd{x 2 P 2 )u{x 2 P 2 ) lv {r 2 -tii Pi)lM+ i- ~ 1)75] 



(P 2 -6Pi) 2 (£iPi-x 2 P 2 ) 2 ' 
define the lowest-order perturbative PS (PT) distribution amplitude and the corresponding 
hard amplitude, respectively. The meaning of the variable £1, regarded as a momentum 
fraction modified by collinear gluon exchanges, will become clear below. 

Consider the 0(a s ) radiative corrections to Fig. 1(a) in the covariant (Feynman) gauge, 
which are shown in Fig. 2, and identify their infrared divergences. Self-energy corrections to 



the internal lines, leading to a next-to-leading-order hard amplitude, are not included. Here 
we summarize only the results of the 0(a s ) factorization, and leave the details to Appendix 
A. It will be shown that all the diagrams in Fig. 2 can be written as the convolution of 
the lowest-order hard amplitudes H^j, in Eq. (jllj) with the 0(a s ) divergent distribution 
amplitudes <j>g T in the collinear region with the loop momentum I parallel to P\. The 
expressions of <pg T will provide a basis of constructing the all-order definitions of the two- 
parton twist-3 distribution amplitudes as nonlocal matrix elements. 

Figures 2(a)-2(c) are the two-particle reducible diagrams with the additional gluon at- 
taching the two valence quarks of the initial state. It has been known that soft divergences 
cancel among these diagrams. The reason for this cancellation is that soft gluons, being 
huge in space-time, do not resolve the color structure of the pion. Collinear divergences 
in Figs. 2(a)-2(c) do not cancel, since the loop momentum flows into the internal lines in 
Fig. 2(b), but not in Figs. 2(a) and 2(c). Inserting the Fierz identity into Figs. 2(a)-2(c), 
we obtain the approximate loop integrals in the collinear region, 

JW,W,(0 W £ /^tf^ne^^l)^^,^), (12) 
n=S,T J 

respectively. The 0(a s ) PS pieces, 

*2<*.«o - f/^*^^*^' -*■> • (13) 

= & I S***' 7 '^^^*'*^ 6 - Xl) ' (15) 

and the PT pieces with 75 in the above expressions being replaced by 75 (yl + fl- — 1), contain 
the collinear (logarithmic) divergences in Figs. 2(a), 2(b), and 2(c), respectively. Note that 
the momentum fraction x\ in Fig. 2(b) has been modified into £1 = x± — 1 + / , because the 
loop momentum I flows through the hard gluon. 

The factorization of the loop integrals associated with the two-particle irreducible dia- 
grams in Fig. 2(d)-2(f) is written as 

I meW) « E f d ^nlne,nM^lM°\^X 2 ), (16) 
n=5,T J 
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with the 0(a s ) PS pieces 

Sd{Xh ^ ~ 2^cv J (2^ d(xijPl)75 (^p7T7F 7 u{xiPl) p—i 

i 



x 



(17) 

- sSe/ ^ lP ^^^^ P 4^-^ ' (18) 

A^ c = 3 being the number of colors. The corresponding PT pieces are defined similarly with 
7 5 in the above expressions being replaced by 75 (^+ fl- — l). Figure 2(g) does not contain a 
collinear divergence. Note that Fig. 2(d) is free of a soft divergence, because the additional 
gluon attaches the hard gluon. The soft divergences cancel between Figs. 2(e) and 2(f). 
The above absence of the soft divergences is obvious from the cancellation in Eq. (|17|) and 
between Eqs. (j!8)) and (fTTIJ) as / — > 0. The contribution from Fig. 2(d) has been split into two 
terms as a consequence of the Ward identity Q . The first and second ^-functions in Eq. (117)) 
correspond to the cases without and with the loop momentum I flowing through the internal 
lines, respectively. The Feynman rule n_^/n_ • I, coming from the eikonal approximation 
(see Appendix A), could be generated by a Wilson line in the direction of n_. This factor 
will not appear, if the proof is performed in the axial gauge n_ ■ A = 0. 

Note that Eqs. (jl7j) - (jl9j) possess different color factors due to different color flows in 
Figs. 2(d)-2(f). Combining the contributions from Figs. 2(d)-2(f), we arrive at 

(/) , 

E r « E / d ^€l^m£\^) , (20) 

i=(d) n=S,T J 

where the PS piece, 



4m J (27r) 4 (xiPi + 1) 1 1 



(21) 



w 

is associated with the collinear gluon emitted from the u quark. The color factor Cp implies 
the factorization of the distribution amplitude from other parts of the process in color space, 
which can be achieved only by summing a set of diagrams. The first and second terms in 
Eq. (pTj) correspond to Figs. 3(a) and 3(b), respectively, where the double lines represent 
the Wilson lines mentioned above. 
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The analysis of Figs. 2(h)-2(k) is similar, and the conclusion is that Fig. 2(h) is split 
into two terms as in Eq. (|17jh the collinear gluons in Figs. 2(i) and 2(j) are eikonalized 
as in Eqs. ()18|) and (jTUj) . respectively, and Fig. 2(k), like Fig. 2(g), does not contribute a 
collinear divergence. The soft divergences cancel among the above diagrams. Summing over 
Figs. 2(h)-2(j), we derive the correct color factor: 



(i) 

i=(h) 



C T J 



X2 



(22) 



n=S,T 



where the PS piece, 



4mp 



¥1°)7 7775 ^75«(xiPi 



-P 



(2vr 



x 



<*(fi - x x ) - 5 ( & - xi + r)+ 



(xiPi - 2 

/+ 



/ 2 n_ • / 



(23) 



is associated with the collinear gluon emitted from the d quark. The first and second terms 
in Eq. (}2*3"j) correspond to Figs. 3(c) and 3(d), respectively. 
The sum of Eqs. O, (JUJ) and (J22J) leads to 



(*) 

i=(o) n=5,T 



(24) 



where <f>^ and are represented by the 0(a s ) terms of the nonlocal matrix elements with 
the structures 75 and 7s(^+ yi- — 1) sandwiched, 

dy 



m 
P+ 



2tt 
dy~ 



m J 2ir 
xPexp 



y 



(0\d(y~) 75 V exp 



'0 



u(0)|u(xP)d(xP)) 



(o|%-) 75 (^ + ^L-l) 



ig J dzn_ ■ A(zu-) 




u{0)\u{xP)d{xP)) 



(25) 



respectively. By expanding the quark field d(y~) and the path-ordered exponential (Wilson 
line) into powers of y~ , the above matrix elements can be expressed as a series of the 
covariant derivatives (D + ) n d(0), and are gauge invariant. The integral over z in fact contains 
two pieces: for the upper Wilson line in Fig. 3(a), z runs from to 00. For the lower Wilson 
line in Fig. 3(b), z runs from 00 back to y~. The light-cone coordinate y~ 7^ corresponds 
to the fact that the collinear divergences in Fig. 2 do not cancel. Notice the different 
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kets \ti + (P)) and \u(xP)d(xP)) in Eqs. (jSJ) and (j2H), respectively. Equation (|23j) plays 
the role of an infrared regulator for parton-level diagrams. A hard amplitude, obtained by 
subtracting Eq. (|25|) from the parton-level diagrams, then corresponds to the regularized 
parton-level diagrams. After determining the gauge-invariant infrared-finite hard amplitude 
H(x), we convolute it with the physical two-parton twist-3 pion distribution amplitudes 
defined in Eq. (jSJ) and graphically shown in Fig. 3(e). Models for the two-parton twist-3 
pion distribution amplitudes have been derived using QCD sum rules |9(. 

It is easy to confirm that (j)$ T is reproduced by the perturbative expansion of the matrix 
elements in Eq. (125)1 . Take Eq. (|21)) as an example. Fourier transforming the gauge field 

A(znJ) into A(l), we have 

f°° ~ n a 

— ig I dzexp[iz(ri- • I + ze)]n_ • A(l) = g — ^—A a {l) . (26) 
Jo n - ' l 

The field A a (l), contracted with another gauge field associated with the u quark, gives the 

gluon propagator in Fig. 3(a). It is then realized that the eikonal propagator is generated 

by the path-ordered exponential. The second piece of the Wilson line corresponds to the 

second term in Eq. ()21j) : 

rv~ n a _ 
— ig I dzexp[iz(n_ ■ I + ze)]n_ • A(l) = —g exp(il + y~)A a (l) , (27) 

J oo n — ' ' 

where the extra Fourier factor exp(il + y~) leads to the function 5(£ — x + l + /Pi )■ The 
field A a (l), contracted with another gauge field associated with the u quark, gives the gluon 
propagator in Fig. 3(b). The Feynman rules for Eq. ()2*3*|) can be reproduced in a similar way, 
where A a (l) is contracted with another gauge field associated with the d quark. Equations 
(fl3 ^1 -(fI3 ]l are derived by contracting the gluon fields associated with the u and d quarks. 

The above derivation also applies straightforwardly to the factorization of the collinear 
divergences associated with the final state, which arise from the region with the loop mo- 
mentum parallel to P 2 . Hence, we have an expression similar to Eq. (|21)l. ^2 m=ST Hm ®0m , 
where is the final-state distribution amplitude, and Hm is the hard amplitude with the 
Fierz identity inserted into the final-state side of H^. The symbol £g> represents the con- 
volution in the variable £2- The momentum fraction associated with the final state will be 
modified into £2 = ^2 — l~ /P2 by ^ ne collinear gluons parallel to P2, if the loop momentum 
I flows through the hard amplitude Hm ■ 

The 0(a s ) radiative corrections to Fig. 1(b) are displayed in Fig. 4. The factorization 
of the collinear divergences from these diagrams is referred to Appendix B. The result is 
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similar to Eq. but without the PT contributions, because of r y u (/i+ fi- — 1)7„ = 0, 
where the gamma matrices 7^ and 7^ come from the gluon vertices in Fig. 1(b). Hence, we 
derive, from Fig. 4, 

(*) , 

|3 I* w y d^^^i^O^fe,^) , (28) 

i=(o) 

where the definition of the 0(a s ) PS distribution amplitude <j)g is the same as in Eq. (|2~K|) . 
This is expected due to the universality. 

In summary, the 0(a s ) factorization of the process rcy* — > tt is written as 

where denotes the complete set of the 0(a s ) corrections, and Hn)n receive the contri- 
butions from both Figs. 1(a) and 1(b) now. The first term on the right-hand side of the 
above expression does not contain the collinear divergences from the loop momentum I par- 
allel to P2. In this region l + is negligible, £1 approaches x%, and the corresponding collinear 
divergences cancel in 0„ 1 . For the similar reason, the second term on the right-hand side 
of Eq. (|29j) does not contain the collinear divergences from I parallel to Pi. That is, the 
initial-state and final-state collinear divergences in have been completely factorized 
into the first and second terms on the right-hand side of Eq. (|29|). respectively. The 0(a s ) 
hard amplitude defined via Eq. (J2HJ), is infrared finite. Note that H^ 1 ' contains the 

contributions from Figs. 2(g), 2(k) and the self-energy corrections to the internal lines. 

Summing Eq. (}2~9*j) and the lowest-order diagrams G^°\ the factorization formula for the 
two-parton twist-3 contributions to the process ivy* — ► rr is given, up to 0(a s ), by 



G (o) + G (D = J- i^ + W) {*$}+<!>$), (30) 

n,m=S,T 

where the trivial factorizations, 

#f=£#S®«, ^=J3^ } ®^. H(1) = E ^®H^^m 

m=S,T n=S,T n,m=S,T 

have been adopted. The last formula in Eq. ()31|) defines the 0(a s ) hard amplitude ifnm-The 
explicit expression of, for example, Hf§, is given by 

2 r 2 / ^[7^757,(A-6 fih^Tsj , tr[Yl 5 l»(Pi-& A)7,7s] ] , 00 , 
H g8 {ti,M - ~ 2 eg ^o| (p2 _ eiPi)2(eiPi _ 6P2)2 + (Pi _ 6 P 2)2( £ lPl _ £ 2 P 2 ) 2 j ^ 32 ) 
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It is obvious that the PS and PT structures must be included simultaneously for a complete 
two-parton twist-3 collinear factorization. 

III. ALL-ORDER FACTORIZATION OF vr 7 * -» vr 

In this section we present the all-order proof of two-parton twist-3 factorization theorem 
for the process 717* — > ir, and construct the gauge-invariant distribution amplitudes in 
Eq. (|25|). It has been mentioned in the Introduction that factorizations of a QCD process 
in momentum, spin, and color spaces require summation of many diagrams, especially at 
higher orders. The diagram summation can be handled in an elegant way by employing the 
Ward identity, 

1»G»{IMM,--- A) = 0, (33) 

where represents a physical amplitude with an external gluon carrying the momentum I 
and with n external quarks carrying the momenta k±, hi, • • • , k n . All these external particles 
are on mass shell. The Ward identity can be easily derived by means of the Becchi-Rouet- 
Stora (BRS) transformation [3^. We shall employ the similar Ward identity, 

llG^{hMMM>--- ,K)i v 2 = o, (34) 

where represents a physical amplitude with two external gluons carrying the momenta 
li and I2, and with n external quarks. 

We shall prove two-parton twist-3 factorization theorem for the process 7r 7 * — > 7r to all 
orders by induction. The factorization of the 0(a s ) collinear divergences has been worked 
out in Sec. II. Assume that factorization theorem holds up to 0(a^): 

G{k) = E EE^®^ _i) ®^> k = 0A,..-N. (35) 

n',m'=S,T i=0 j=0 

Q( k ) denotes the parton-level diagrams of 0{a k s ) with G^ shown in Fig. 1. The initial-state 
distribution amplitude <j>®, is defined by the 0{a\) terms in the perturbative expansion of 
Eq. (J25)) . and the final-state distribution amplitude (p^, defined similarly by the complex 
conjugate of Eq. (|25jl. H^^T is the remaining 0(a;^ _l_J ) piece of the process, which does 
not contain the infrared divergences. Equation (|35|) implies that all the initial-state and 
final-state collinear divergences in G^ have been collected into (jy^j and 4>~$ systematically. 

13 



Inserting the Fierz identity, we also obtain the trivial factorizations of the distribution 
amplitudes </> and the diagrams G at arbitrary orders of a s , similar to Eq. (J3*T|) . We then 
have the factorization, 

g&= E EE^l>®H^r j) ®^ m , (36) 

n',m'=S,T i=0 j=0 

in which, for example, <j)gg contains the piece (fissu extracted from Eq. (J5TJ, 



w*i,w - 4mo y (27r ^ 



^1 A+ / ^ 

{XiPi + I) 2 



1 n_ 



73 



P n_ • / 



5(6 - a?i) - ^ ( Cl - a?i + 7TT 



^ + yj ■ (37) 

Below we prove the collinear factorization of the 0(a^ +l ) diagrams 

G (N+1) 

assuming 

Eq. ()35jl or (|36p. Look for the radiative gluon in a subset of 0(af +1 ) diagrams G^ N+1 \ one 
of whose ends attaches the outer most vertex on the upper u quark line. Let a denote the 
outer most vertex, and j3 denote the attachments of the other end of the identified gluon 
inside the rest of the diagrams. There are two types of collinear configurations associated 
with this gluon, depending on whether the vertex (3 is located on an internal line with the 
momentum along P\. The fermion propagator adjacent to the vertex a is proportional to 
Pi in the collinear region with the loop momentum I parallel to P\. If /3 is not located on a 
collinear line along Pi, the component 7 + in 7° and the minus component of the vertex (3 
give the leading contribution. If (3 is located on a collinear line along Pi, f3 can not be minus, 
and both a and (3 represent the transverse components. This configuration is the same as 
of the self-energy correction to an on-shell particle. According to the above classification, 
we decompose the tensor g a p appearing in the propagator of the identified gluon into 

9a/3 = S a+ Sf3- - 5 a _L5p_L + 5 a -5f3 + . (38) 

The first (second) term on the right-hand side of Eq. (|3*K|) extracts the first (second) type of 
initial-state collinear divergences mentioned above. The third term does not contribute due 
to the equations of motion in Eq. (JBJ) . 

We discuss the factorization of the first type of collinear configurations denoted by G\ N+1 ^ . 
As stated before, the identified collinear gluon with a = + and (3 = — does not attach the 
upper or lower quark line directly, which carries the momentum along P\. That is, those 
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diagrams with Figs. 2(a) and 2(b) as the 0(a s ) subdiagrams are excluded from the set of 
q( n + 1 ) as discussing the first type of collinear configurations. We employ the replacement, 

<W V - , (39) 

where the light-like vector n_ a also selects the plus component of 7°, and l@ selects the 
minus component of the vertex j3 in the collinear region. lp can attach all the internal lines, 
no matter they are or are not parallel to Pi. When it attaches a line parallel to Pi, the result 
diminishes. That is, Eq. indeed picks up the first type of collinear configurations. We 
then consider Fig. 5(a), which contains a complete set of contractions of lp, since the second 
and third diagrams have been included. The contractions of lp, represented by arrows, hint 
the application of the Ward identity in Eq. (JHHj) to the case, in which the on-shell external 
u quark, d quark and gluon carry the momenta £iPi, x\P\ and I, respectively. The Ward 
identity states that the expression in Fig. 5(a) vanishes. The second and third diagrams in 
Fig. 5(a) give 

h ? p y A(eifi) = J I , (/-& A + 6 hXllPl) = -U{1 X P X ) , 
U A" / « A" / 

IpdixxPx)^ 1 = -d{x x P x ) , (40) 

respectively, where the equations of motion in Eq. (JHJ) have been employed. The terms 
ii(£iPi) and d{x\P\) at the ends of the above expressions are associated with the 0(a^) 
diagrams. 

We then insert the Fierz identity into Fig. 5(a), and factor the lowest-order expressions 
(i(xiPi)rM(^iPi) with T being the PS or PT structure considered in this work. The result is a 
relation shown in Fig. 5(b), where the cuts on the quark lines denote the insertion of the Fierz 
identity, and the double (Wilson) lines represent n_ Q /n_ ■ I in Eq. (jHOJ)- Figure 5(b) implies 
that the diagrams G\ N+1 ^ corresponding to the first term in Eq. (JHHj) are factorized into the 
convolution of the full diagrams with <pnu, n — S,T. The same discussion applies to 
the factorization of the diagrams Gj Ar+1 ' ) with the collinear gluon emitted from the outer 
most vertex on the d quark line, leading to a convolution of Gn with (jy'j. Similarly, for 
the subset of G^ N+1 \ in which the identified radiative gluon emitted by the outgoing quarks, 
the collinear divergences are also classified into two types. In this case it is the third term 
in Eq. (jHHj) that corresponds to the first type of collinear divergences, and the replacement 
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in Eq. is modified into 



6 a -6 p+ - . (41) 



The pieces cf)mu and containing the identified gluon are factored out of G\ N+1 ^ in the 
collinear region with the loop momentum parallel to P^. 

We conclude that G^ +1 \ in which the replacement in Eq. ()39|) is applicable to the initial- 
state side, and , in which the replacement in Eq. (|4*T]) is applicable to the final-state 
side, are written as, 



n=S,T 

- E c^Sr («) 



m=5,T 

respectively, where ^Lau represents 

0« ... = 0« + . (44) 

^n(m)|| ^n(m)u 1 ^n(m)d V / 

Equation (j4*2*)l is displayed in Fig. 6. Other diagrams, which do not contain the radiative 
gluons on either the initial-state side or the final-state side, are self-energy corrections to 
internal lines. They are infrared finite, and contribute to the 0(a^ +1 ) hard amplitude. 

The above procedures are also applicable to the 0(a{ +1 ) initial-state and final-state 
distribution amplitudes <j>n and <pm ■ We identify the gluon in the complete set of 
0(ai +1 ) diagrams <fin , one of whose ends attaches the outer most vertex a on the u quark 
line. The other end attaches the vertex (3 inside the rest of the diagrams. For the first term 
on the right-hand side of Eq. (J38|) . we have a Ward identity similar to Fig. 5(a). Figure 5(b) 
then leads to the factorizations of the initial-state and final-state distribution amplitudes, 



n\\ ~ 

n'=S,T 

m\\ ~ 

m'=S,T 



E^W^- (45) 

i'=S,T 

E ® €w , (46) 

n'=S,T 

where the PS and PT structures in Eq. Q have been inserted. 

We sum Eqs. (142)1 and ()43)1 . and subtract the double-counted diagrams G\^ +1 \ to which 
the replacements are applicable to both the initial-state and final-state sides as shown in 
Fig. 7(a). Note that Gy +1 ' does not contain the diagram in Fig. 7(b), in which the same 
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n,m=S,T i=0 j=0 



gluon is identified in G^ +1 ^ and simultaneously. This type of diagrams are not 

double-counted, and Eq. (JSjj) does not apply to them. The factorization of the first type of 
collinear divergences associated with the identified radiative gluon from G\ N+1 \ N > 1, is 
then given by 

MN+l) n (N+l) MN+1) MN+1) 
I — II" II — II II ' 

- E Kli ® G $z ® ^ + ^ ® G ^ ® <i - *S ® ® ^ii] » ( 4T ) 

n,m=S,T 

where the trivial factorizations, similar to Eq. have been inserted. For the factorization 
of we rely on the Ward identity in Eq. (JH3). Substituting Eqs. flHHD, (jEJ), and flUJ) 

into Eq. (|4"7j). a simple algebra gives 

AT iV-i 

G f +1) - E E E [«C ® H ^~ j) ® *SP + ^ ® ® «C 

(48) 

It is obvious that Eq. (J48j) is not Lorentz covariant, since G\ N+1 \ <pty and include only 
the first or third term in Eq. ((HHJ)- The Lorentz covariance can be recovered by adding the 
other two terms in Eq. (|H8|) into the tensors for the identified radiative gluons on both sides 
of Eq. (|48|). As recovering the Lorentz covariance, the second type of collinear configurations 
associated with the tensor — 5 a ±5p± is taken into account. We then further add the infrared- 
finite (self-energy correction) diagrams mentioned above. At last, the complete collinear 
factorization of the diagrams 

G (N+l) 

is given by 

N N-i 

g {n+1) = E E E M m) ® H ^~ j) ® *$? + $ ® ® <e x) 

n,m=S,T i=0 j=0 

® ^- J '- 13 ® <^ +1) ] + , (49) 

where the 0(a^ +1 ) function F^ N+1 ^ contains the contribution corresponding to the difference 
between 8 a +8p- and n_ a lp/n- ■ I, except the infrared-finite diagrams. Equation (|49|) can be 
simplified into 

N+l N+l-i 

GiN+ i } = j2 J2 E ® ® « , (so) 

n,m=S,T i=Q j=0 

with the 0(a^ +1 ) hard amplitude H^m~^ being defined via 

F (* + D = J2 ® ff£ +1 > ® . (51) 



n,m=S,T 
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Equation (jSUJ) states that all the two-parton twist-3 collinear divergences in the process 
7T7* — > 7r have been factorized into the distribution amplitudes in Eq. (J2"5|) order by order. 

Before closing this section, we prove that soft divergences do not exist in the process 
7T7* —>■ 7i at the two-parton twist-3 level. The 0(a s ) soft cancellation has been demonstrated 
in Sec. II. Assume that the O(a^) diagrams G^ N ' do not contain any soft divergences, 
though they contain collinear ones. They are then dominated by hard dynamics, by collinear 
dynamics associated with the initial-state pion, and by collinear dynamics associated with 
the final-state pion. Consider the 0(a^ +1 ) diagrams G^ N+1 \ Similarly, we look for the 
radiative gluon radiated from the outer most vertex on the u quark line in the initial state, 
and adopt the decomposition of the tensor g a/3 in Eq. The attachment of a soft gluon 

to an off-shell internal line does not introduce an infrared divergence, since an off-shell 
propagator behaves at least like 1/Q. As the soft gluon attaches a collinear line along 
P2, the vertex (5 must be dominated by the minus component. The gamma matrix 7" is 
dominated by the component 7" 1 " stated above. Therefore, the replacement in Eq. ()39)1 still 
holds for the first term on the right-hand side of Eq. (|3*K|) . Similarly, the second term on the 
right-hand side of Eq. (|38jl corresponds to the attachment of the soft gluon to a collinear 
line along Pi. Again, the third term on the right-hand side of Eq. (}3*Hj) does not contribute 
due to the equations of motion. 

The above reasoning indicates that the configurations associated with the identified soft 
gluon are the same as those associated with the identified collinear gluon. The procedure for 
deriving the collinear factorization applies, and the soft divergences in the process Try* — > tt, 
if there are any, can be absorbed into the initial-state and final-state distribution amplitudes. 
According to Eq. (|33j). the distribution amplitudes (pi for j = 1, 2, • • • , N are free of soft di- 
vergences by assumption. We have Eq. Po"]) for the factorization of the 0(a^ +1 ) distribution 
amplitudes 4>^ +1 ^ ■ Because the soft divergences cancel among the diagrams for <fity shown in 
Fig. 3 [between Figs. 3(a) and 3(b) and between Figs. 3(c) and 3(d)], and do not contain 
soft divergences, <f)^ +1 ^ do not either. We then turn to the 0(a^ +l ) distribution amplitude 
associated with the tensor —5 a ±5p±. If it contains soft divergences, the evaluation of these 
divergences, arising from an incomplete set of diagrams (the identified soft gluon does not 
attach the Wilson line), is not gauge invariant. This result contradicts the gauge-invariant 
definition of the distribution amplitudes in Eq. ()25|) . Hence, the 0(a^ +1 ) distribution am- 
plitude must be free of soft divergences. Extending N to infinity, the absence of the soft 
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divergences in the distribution amplitudes 0„ jm and in the diagrams G is proved. We then 
complete the all-order proof of the two-parton twist-3 collinear factorization theorem for the 
process 717* — > ir. 



IV. FACTORIZATION OF B — ► ttIu 



As emphasized in the Introduction, the contributions to the exclusive semileptonic decay 
B — > -nlv from the two-parton twist-3 pion distribution amplitudes <ps,T are of the same 
power as those from the twist-2 one. These contributions are originally proportional to the 
ratio itlq/Mb- However, the corresponding convolution integral for the B — > 7r form factor is 
linearly divergent in collinear factorization theorem, such that it is proportional to the ratio 
Mb/A, if regulated by an effective cut-off x c ~ A/ Mb- Combining the two ratios mo /Mb 
and Mb/A, the two-parton twist-3 contributions are in fact not down by a power of 1/Mb'- 

— / dx 2 H s {x 2 )<ps{x 2 ) ~ rrz- I dx 2 — ~ O (-=-) , (52) 

for the asymptotic functional form (f>s ~ 1 [9] and if^ ~ 1/^1 m ^Q- ^3|) below. The 
presence of the linear divergences modifies the naive power counting rules, and the two- 
parton twist-3 contributions become leading-power in the B meson transition form factors 
~l 

24| . This is our motivation to prove the corresponding factorization theorem. Note that 
the two-parton twist-3 contributions are next-to-leading-power in the pion form factor. 

The momentum Pi of the B meson and the momentum P 2 of the outgoing pion are 
parameterized as 

Pi = ^|(l,l,0 T ), P 2 = ^(0 )?7 ,O T ), (53) 

where rj denotes the energy fraction carried by the pion. Consider the kinematic region with 
small q 2 , q = P\ — P 2 being the lepton pair momentum, i.e., with large 77, where PQCD 
is applicable. Let the light spectator d quark in the B meson (pion) carry the momentum 
ki {k 2 = x 2 P 2 ), and the b quark carry the momentum Pi — k\. We have the equations of 
motion, 

(A- h ~ m b )b{P 1 -k x ) = 0, d(h) ft = , (54) 
where the d quark mass ro^ has been neglected, and those for the valence quarks in the pion. 
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A. 0{a s ) Factorization 



We start with the 0(a s ) collinear factorization of the final-state distribution amplitudes. 
The lowest-order diagrams for the B — > irW decay are the same as in Fig. 1, but with 
the upper quark line in the initial state representing a b quark and with the symbol x 
representing a weak decay vertex. Figure 1(a) gives the amplitude, 

„(o), \ -g 2 c F u(x 2 p 2 )Y(p2- h)lA p i - k 1 )d(k 1 hJ(x 2 P2) 

G ^^)-— (p 2 -k 1 nk 1 - X2 p 2 r ' (55) 

where the higher-power term ]k\ in the numerator has been dropped. Inserting the Fierz 
identity in Eq. @ into the above expression, we obtain the trivial factorization formula, 

G (0) (xi,x 2 ) = J ^4 0) (zi,6)4 0) (z2,6) , (56) 

with the lowest-order hard amplitude and distribution amplitude, 

Hy (xi,&) = — ~ — m - 



13 " ~^ m ° (p 2 - mh - X2 p 2 + iy 



-g 2 C F tr[Y ~ kiMhh^ 5 } 

-m - 



x 2 ,i 2 ) = -^—u{x 2 P 2 )-i b d(x 2 P 2 )5(i 2 -x 2 ). (57) 



5,(0) < 1 
s 4m 



The PT structure does not contribute because of •y lJ (f l + A- ~ tylslu = 0. It is observed 
that depends only on the single plus component of k\ through k\- P 2 , which defines the 
momentum fraction x\ = k^ / P^ carried by the d quark in the B meson. 

Consider 0(a s ) corrections to Fig. 1(a), which are displayed in Fig. 4 with the initial and 
final states being flipped. Here we summarize only the results of their factorization, and 
refer the details to Appendix C. The factorization of the two-particle reducible diagrams 
in Figs. 4(a)-4(c) is straightforward. After inserting the Fierz identity, the loop integrals 
associated with Figs. 4(a), 4(b), and 4(c) are written as 

I iaUbUc) « / ^2Hf\x u ^ 2 )4l SbtSc (x 2 ^ 2 ) , (58) 
where the PS collinear pieces are given by 

A 1 )/ c\ ~ i 9 2G F f dH x 2 P 2 + / p 1 -j 4 , , , . 

(1) i 9 2 c F r dH x 2 p 2 + / x 2 p 2 -/ - f p A,f, r\ 

^ (X2 ' 6)= 1^77 J^y u{X2P2)l \x 2 p 2 + iy l5 {x 2 p 2 - i)^ d ^p 5 ^ 2 ~ X2 + pff } 
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The momentum fraction £2 in Fig. 4(b) has been modified into £2 = x 2 ~ l~ /P 2 by the 
collinear gluon exchange. 

The collinear factorization of Figs. 4(d)-4(g) is summarized as 

I (e) + I {9) ^ f<% 2 Hf ( Xl , &)4l,{x 3 , 6), (62) 

with the PS pieces, 

,(i)/ c n -W 2 f dH (- d n a x 2 P2+ f J/ px 1 

</>sJ (^2,6 =0 — 7T" / 7r-u u a; »^ T F p T n2 ^d{x 2 P 2 

2m C F J (27r) 4 (^2^2 + t) l 2 n + -l 



x 



5(6 - ^2) - 5 ( & - x 2 + — 



P-2 



(63) 

+^>ton&J ^^^^^h^i 6 ^ • (64) 

The contribution from Fig. 4(d) has been split into two terms as a consequence of the Ward 
identity. The first and second terms correspond to the cases without and with the loop 
momentum / flowing through the hard gluon, respectively. The Feynman rule n + p/n + ■ I in 
the collinear divergent pieces, coming from the eikonal approximation, can be represented 
by a Wilson line in the direction n + . Note that each of the hard amplitudes from Figs. 4(e) 
and 4(g) contains a residual dependence on I. It is their sum that does not depend on I, and 
appears in the desired factorization form shown in Eq. (|62[) (see Appendix C). 

Equations ([63 )1 - ([65 )1 carry different color factors due to different color flows in Figs. 4(d)- 
4(g). Their sum leads to the factorization in the color space: 

(9) , 

I' ~ / dt 2 Hf\x x ^ 2 )cj>f u {x 2 ^ 2 ) , (66) 

i=(d) 

where the PS collinear piece, 

rhWfr P \ - ~ 1 ^ Cf f dH .,(- P ^ Il^LL. Mt P ^ ™ +/3 
9Su\ X ^Kl) - ~ A / /o xa U K X ^2)1 1- p I A2 75»l^2-f2 



4m J (2tt) 4 y " (x 2 P 2 + l) 2 ' y J Pn + -l 
<K& -x 2 ) -5 (^ 2 - x 2 + 



(67) 



is associated with the collinear gluon emitted from the u quark. The appropriate color factor 
Cf indicates the factorization between the distribution amplitudes and the hard amplitude 
in color space. 
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The collinear factorization of Figs. 4(h)-4(k), derived in a similar way, is written as 

(k) 



f^P ~ [ d^ 2 Hf{x u 6)0^(^2, 6) , (68) 



where the PS piece, 



X2l6) " 1^7 J (2^ M(X2 2)75 (x 2 p 2 - if 1 d[X2P2) 



4m J (2tt) 4 v z ZJIO (x 2 P 2 -iy 1 " z ZJ Pn + -l 

(69) 



x 



<K6 - z 2 ) - * ( 6 - ^2 + T^r 



is associated with the collinear gluon emitted from the d quark. Note that Fig. 4(k) also 
contributes a collinear divergence, and we have to combine the results from Figs. 4(j) and 
4(k), so that the hard amplitude does not depend on the loop momentum I. At last, the 
sum of Eqs. ()58|) . (JB5J) and ((QHl) gives 

(*) , 

J2 « y dC 2 Hf\ Xl , &)#(a; 2 , 6) , (70) 

i=(o) 

where the PS collinear piece </>^ is defined by the 0(a s ) term of the complex conjugate of 
Eq. (|25|). consistent with the universality. 

The amplitude corresponding to Fig. 1(b) is written as 

(0) -g 2 C F u(x 2 P 2 )-f^(P 1 -x 2 P 2 + m b )YK p i - k 1 )d(k 1 )-f u d(x 2 P 2 ) 

lXl ' X2j " 2 [(P x - * 2 P 2 )s - mgPi - * 2 P 2 ) 2 ' l j 

which, after inserting the Fierz identity, leads to the trivial factorization formula, 

Gi ° } = E [t&fflfri, 6)^(^2, 6). (72) 
The lowest-order hard amplitudes and distribution amplitudes of the PS and PT structures 



are written as, 



g 2 C F tr[ lfi ( A - 6 A + m h )rKPi - fci)d(fcih„75] 



[{Pi ~ £2 A) 2 - mg] (Ar x - 6 A) 



-<? 2 Cf M7m( A - 6 A + m b )Yb(Pi - ^1)^(^)7,75] 



2 



r (o), t s . -<? 2 CV tr[ 7M ( A - 6 A + m b )Yb(Pi - fc 1 )d(fc 1 ) 7l/ (^+ yL - 1)75] 



#T Vl>6) = o m 



2 u [(A-kA^-mgK^-kA) 2 



<.(«)/ 



1 



-r (^2,6) = ~ A u(x 2 P 2 )>y 5 (fl + ft- - l)d(x 2 P 2 )5(£, 2 - x 2 ) . (73) 

4mo 
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The lowest-order PS distribution amplitude ( X2 ^ 2 ) is the same as in Eq. (|57J). 

Below we discuss the collinear divergences in the 0(a s ) corrections to Fig. 1(b) in the 
covariant gauge, which are displayed in Fig. 2 with the initial states and the final states 
being flipped. The details are referred to Appendix C. Figures 2(a)-2(c) are factorized 
straightforwardly, leading to 

j(.),W,(c) „ J2 f d^\x u ^ a>mbtmc (x 2 ^ 2 ) . (74) 

m=S,T J 

The PS collinear divergent functions <^ Sh Sc (x 2 , £2) are the same as those shown in Eqs. (fHHJl 
- respectively. The PT pieces have the similar expressions with 75 being replaced by 

75 {i+ i- - 1)- 

For the irreducible diagrams Figs. 2(d)-2(g), a summation of their contributions is nec- 
essary for obtaining the desired collinear factorization. Figure 2(d) is split into two terms 
after applying the Ward identity, similar to Eq. (JHBJ). The sum of the first term and Fig. 2(e) 
gives the correct color factor, which corresponds to the separate color flows between the dis- 
tribution amplitudes and the hard amplitude. However, the hard amplitude still depends on 
the loop momentum I, which is not yet in the desired form. We must further add Fig. 2(g) 
to this sum in order to obtain the hard amplitude without the I dependence. The sum of the 
second term from Fig. 2(d) and Fig. 2(f) also gives the correct color factor, and the result 
represents the case with / flowing through the hard amplitude. Hence, we arrive at, from 
the combination of Figs. 2(d)-2(g), 

(9) . 

Yl d ^K^ GMSlfa, 6). (75) 

i=(d) m=S,T J 

The collinear factorization of Fig. 2(h)-2(k), derived in a similar way, is written as 

(*) , 

E /^2Fi 0) (^1,6)0^2, 6), (76) 

i=(h) m=S,T J 

with the PS pieces <j}g\ and 4^g\ shown in Eqs. ()67j) and ()69|) . respectively. At last, the sum 
of Eqs. (f74]l . (|75 jl and (|76jl gives 



J2 P * E dbH^xi&WSfab) , (77) 

i=(a) m=S,T J 
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where the collinear divergent functions 4>m\x2,£,2) are defined by the 0(a s ) terms of the 
complex conjugate of Eq. (|23j). 

The factorization of the soft divergences from the decay B — > -nlv has been per- 



u, 



formed in |1|, which results in two B meson distribution amplitudes and <f)^ 



HQ 



arising from the insertion of the fourth and fifth terms of the Fierz identity on the initial-state 
side: 



hjhk = -^hkhj + -^{la)ik{l a )ij + 2 



ik 



1 

+ 4 



IT! 75 fi+ ^ + 1 " 



V2 



75 M/' + l) 



i/c 



v/2 



(78) 



The Wilson line on the light cone can be constructed, only if the hard scale for exclusive B 
meson decays is of O ( a/AM#) . This Wilson line is crucial for the gauge-invariant definitions 
of the distribution amplitudes as nonlocal matrix elements. 

Following the similar procedures in Sec. II, we derive the 0(a s ) factorization of the decay 
B — > ttW, 



(79) 



m=S,T 



where the hard amplitude receives the contributions from both Fig. 1(a) and 1(b). 
Consequently, the factorization formula for the two-parton twist-3 contributions is written, 
up to 0(a s ), as 



G (o) + G (D = £ &$> + )) ® (F£ + flffi 



, ^ (0(0) + Aih 



n=+,- 
m=S,T 



The definitions for the hard amplitudes Hnm and Hnm are similar to those in Eq. 
example, the explicit expression of HYg is given by 



(80) 



D. For 



-g 2 C F j trtfifr- h)ly.{h + M B ) ^_ 7 57,7 5 ] 



-m 



2V2 """{ (P2-k 1 ) 2 (k 1 -x 2 P 2 + l) 2 

trMfi - 6 A + m 6 ) 7 y ( > P 1 + Mg) ^-7 5 7,7 5 ] 
[(Pi-k^-m 2 ]^-^) 2 



(81) 
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B. All-order Factorization 



The all-order proof of the two-parton twist-3 factorization theorem for the process rcy* — > 
7r in Sec. Ill can be generalized to the B — > nW decay with minor modifications. Here we 
highlight only the different points of the proof. In the case of B meson decays there is 
no collinear divergence associated with the initial state, since the b quark is massive, and 
the light spectator d quark is soft [1]. Hence, the B meson side is dominated by the soft 
divergence. The collinear configurations associated with the final-state pion are the same 
as in the process tcj* — > tt discussed in Sec. II. The important infrared divergences are then 
classified into the soft type with small loop momentum I and the collinear type with I parallel 
to P 2 . We shall compare the factorizations of the soft divergences into the initial state and 
of the collinear divergences into the final state. 

Identify the soft gluon emitted from the outer most vertex a on the b quark line in the 
0(af +1 ) diagrams G^ N+1 \ Let (3 denote the attachments of the other end of the identified 
gluon inside the diagrams. The attachment of the soft gluon to collinear lines along P 2 and 
to soft lines gives soft divergences. The soft lines include the b quark line, the spectator d 
quark line, and soft internal lines. The attachment to hard lines off-shell by 0{{P 2 — ki) 2 ) ~ 
0(KMb) also gives soft divergences When the identified gluon attaches the collinear 
lines along P 2 and the hard lines along P 2 — k±, the vertex (3 inside is mainly minus, 
and the vertex a on the b quark line is mainly plus. Therefore, the replacement in Eq. (|39|) 
still works. Certainly, Eq. ()H9j) changes the soft divergences from the attachment of the 
identified gluon to the soft lines. However, these soft lines appear only on the B meson 
side, and will be compensated, when recovering the Lorentz covariance, a case similar to the 
factorization of the second type of collinear configurations in the process 7T7* — > n. 

For the identified soft gluon emitted by the spectator d quark, the first term 8 a+ 8p- on 
the right-hand side of Eq. (}3*H|) corresponds to the configuration, where the identified gluon 
attaches all types of lines except the spectator d quark line. This configuration is the same as 
that associated with the soft gluon emission from the b quark stated above. That is, Eq. ()H9j) 
holds. The second term — 8 a ±8p± corresponds to the configuration with the identified gluon 
attaching the soft lines. The third term 8 a -8p + does not contribute because of the equations 
of motion. 

Since the collinear configurations for the final-state side are the same as those in the 
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process 7T7* — > tt, Eq. ()41|) holds. Based on the above discussion, the symbol || in B — > irlv 
has the same meaning as in 717* — > tt, and the treatment of G\ N+1 ^ follows the procedure in 
Sec. Ill [through Figs. 5(a) and 5(b)]. The above reasoning also applies to the factorization 
of the 0(ai J+1 ' 1 ) initial-state distribution amplitudes 0„i| ■ 111 summary, we obtain the 
similar factorizations but with the subscript n = +, — for the initial state 

Following the same induction procedure, we obtain the factorization of the full diagrams 

G (iY + l) ; 

N+l N+l-i 

G {N+1) = E E E <® ® ® « > (82) 

"= + ,- 1=0 7=0 

m=S,T 

where the 0(a^ +1 ) hard amplitude i?nm +1 ' ) is infrared finite. Equation (|82)l indicates that 
all the soft and collinear divergences in the semileptonic decay B — > irlv can be factorized 
into the distribution amplitudes <pn and 4>m at the parton level order by order, and that 
the proof of the corresponding two-parton twist-3 factorization theorem is completed. These 
parton-level distribution amplitudes serve as the infrared regulators for the derivation of the 
hard amplitudes from the parton-level diagrams. We then compute the B — ► tt transition 
form factor by convoluting the hard amplitudes with the meson distribution amplitudes, 
in which the quark states are replaced by the physical B meson and pion states. Both 
the twist-2 and two-parton twist-3 contributions to the B — > tt form factors F + (q 2 ) and 
F (q 2 ) in the standard definition have been evaluated in 24]. It has been observed that the 
latter are of the same order of magnitude as the former, consistent with the argument that 
the two-parton twist-3 contributions are not power-suppressed and chirally enhanced. The 
light-cone sum rules also give approximately equal weights to the twist-2 and two-parton 
twist-3 contributions to F + j^]. 



V. CONCLUSION 



In this paper we have investigated the infrared divergences in the process 7T7* — >• it at 
the two-parton twist-3 level. We summarize our observations below. There are no soft 
divergences associated with the pion, since they cancel among diagrams. The absence of the 
soft divergences is related to the fact that a soft gluon, being huge in space-time, does not 
resolve the color structure of the color-singlet pion. In the collinear region with the loop 
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momentum parallel to the pion momentum, we have shown that the delicate summation of 
many diagrams leads to the 0(a s ) factorization in the momentum, spin and color spaces. 
We have presented an all-order proof of the two-parton twist-3 factorization theorem for 
the process Tf-f* — > it by means of the Ward identity. This proof can also accommodate the 
twist-2 factorization theorem presented in |l| and the twist-4 one simply by considering the 
corresponding structures in the Fierz transformation in Eq. (JSJ). 

The idea of the proof is to decompose the tensor g a p for the identified collinear gluon 
into the longitudinal and transverse pieces shown in Eq. (|38|). The longitudinal (transverse) 
piece corresponds to the configuration without (with) the attachment of the identified gluon 
to a line along the external momentum. The former configuration can be factorized using 
the Ward identity as hinted by the replacement in Eq. (J39j) . The factorization of the latter 
configuration can be included by demanding the Lorentz covariance of the factorization. 
We emphasize again that the parton-level distribution amplitudes serve as the infrared 
regulators for the derivation of the hard amplitudes from the parton-level diagrams, similar 



to the effective diagrams drawn in SCET [2i 



271 ] . The hard amplitudes are then derived 



from the "matching procedure" . Based on the perturbative construction of the distribution 
amplitudes, we have derived their two-parton twist-3 definitions as nonlocal matrix elements, 
where the path-ordered Wilson line appears as a consequence of the Ward identity. Note 
that our technique is simple compared to that based on the "A-forest" prescription in [5], 
and explicitly gauge invariant compared to that performed in the axial gauge |2j. 

We have generalized the proof to the more complicated semileptonic decay B — > irlv. 
The collinear factorization for the final-state pion is the same as in the process 7T7* — > it. 
The identical collinear structures in both processes justify the concept of universality of 
hadron distribution amplitudes in PQCD. The factorization of the soft divergences for the 

n 

initial-state B meson has been discussed in The conceptual differences are summarized 
as follows. The attachments of a soft gluon to the hard lines off-shell by O(AMg) lead to 
soft divergences. These divergences, like those from the attachments of a collinear gluon 
to the hard lines in the process it^* — > tt, are crucial for constructing the Wilson line, that 
guarantees the gauge invariance of the B meson distribution amplitudes. This explains why 
the characteristic scale of exclusive B meson decays, if factorizable, must be of O(AMg). The 
decomposition of the tensor g a p for the identified soft gluon in Eq. (|3*Hj) and the replacement 
in Eq. (|3*H|) still work. The procedures of the proof then follow those for the pion form factor. 



27 



For a practical application to the B — > ttIv decay, the parton transverse momenta kr mus ; 



be taken into account in order to smear the end-point singularities in the hard amplitudes 



381 ] . This observation implies the necessity of proving k? factorization theorem |T 7, 18]. The 



lplitudes 

HQ. 
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proof of kx factorization theorem is basically the same as proposed in this paper: we simply 
retain the dependence on the loop transverse momenta in hard amplitudes [19| . The relative 
importance of the twist- 2 and two-parton twist-3 contributions to the B — > 7r transition form 
factor has been investigated in [24j, which confirms our motivation to prove the two-parton 
twist-3 factorization theorem: the latter contributions are not power-suppressed and chirally 
enhanced. In a future work the proof will be generalized to nonleptonic B meson decays, such 
as B — > Kit and itit j^Jl^c]]. The corresponding factorization theorem is more complicated, 
since nonleptonic decays involve three characteristic scales: the W boson mass Mw, Mr, 
and the factorization scale of 0(A), such as the parton transverse momenta kr 42 1. 



At last, we compare our construction of collinear factorization theorem in perturbation 
theory with that in SCET. In the former one starts with Feynman diagrams in full QCD. 
Look for the leading region of the loop momentum defined by Eq. (jl)). in which one makes 
the power counting of the Feynman diagrams. It can be found that the approximate loop 
integral in the leading region is represented by the diagram in Fig. 3(e), which leads to the 
definition of a distribution amplitude. In SCET one first constructs the various effective 
degrees of freedom describing infrared dynamics and the effective interactions, and defines 
their powers. Select a specific effective operator, such as those nonlocal operators in Eq. Q. 
Draw the diagrams based on the effective theory, and then make the power counting. It can 
be shown that the diagram in Fig. 3(e) scales like the selected operator, and builds up the 
distribution amplitude. It is not necessary to analyze the infrared divergences in diagrams 
at this stage. That is, one arrives at Fig. 3(e) through approximating loop integrals in the 
full theory in PQCD, but does at the operator and Lagrange level in SCET. Despite of 
the different reasonings for deriving a collinear factorization formula, the calculation of the 
Wilson coefficients is the same. As calculating the Wilson coefficients associated with the 
effective operators from the matching procedure in SCET, the infrared divergences need to 
be analyzed, and their cancellation between the full theory and the effective theory must 
be demonstrated explicitly. This procedure is in fact the same as the derivation of the 
hard amplitudes (Wilion coefficients) in PQCD, where the subtraction of the distribution 
amplitudes (the effective theory) from the parton-level diagrams (the full theory) is done. 
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Therefore, it is legitimate to claim that the constructions of collinear factorization theorem 
are equivalent between PQCD and SCET Q. 
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APPENDIX A: 0(a s ) CORRECTIONS FROM FIG. 2 

In this Appendix we present the details of the collinear factorization of Fig. 2 associated 
with the initial state. The loop integral of Fig. 2(b) is given by 

I (b) = ~eg A C 2 F J ^ i J(xiPi)7 /3 ^^ >1 _ ^ a d{x 2 P 2 )u{x 2 P 2 ) 
Inserting the Fierz identity in Eq. (J5J), we obtain 

m „ w 2 c F f dH - Xl A- / Xl A+ / p ( - P xi 

% 2n tThJ(x 2 p 2 )u(x 2 p 2 ) 1 a (p 2 -x 1 a+ DiM 



2" r u (P 2 -x 1 P 1 + l) 2 (x 2 P 2 -x 1 P 1 + l) 2 
ig 2 C F f dH -j, xi A- / / / / A+ / a /_ px 1 



4m 7 (2tt) 4 v x a/ (anPi - 2 ~ (^iA + 2 ^ 2 

tr[ 7a J(x 2 P 2 ) W (x 2 P 2 ) 7 a (A-^i A+ frUfi+ fi-- 1)75] 
2~* ~ r "" J (P 2 - x a Pi + Z) 2 (x 2 P 2 - ziA + 2 



To derive the above expression, the twist- 2 structure (75 fi-)ih{ fi+lsjji has been dropped. 
The dependence on l~ and on It in and in being subleading according to Eq. (jlj), 
needs to be neglected. Inserting the identity J g?£l<K£i ~~ x i + ^ + /A + ); the first factors of 
the above two terms on the right-hand side of Eq. (jA2|) give the collinear divergent pieces 
"' £1,61) and 0rft(xi,^i) defined in Eq. (JHJ). The collinear factorization of Figs. 2(a) and 



y 56 



2(c) can be performed in a similar way, leading to Eqs. (fTH|) and (|15|) . 
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The loop integral from Fig. 2(d) is written as 

7 = TTTT / T^I^iAjT d{x 2 P 2 )u{x 2 P 2 )^ — - ■ n2 7^- 



ti{T c T b T a )Y cb g a 

^^ PM^yM^ftP ' (A3) 

with the color matrices T a ' b ' c and the triple-gluon vertex, 

r?/? a = / cfca [^/3(2/ - SiPi + x 2 P 2 ) A + ^(2^ - 2a: 2 P 2 - Z) Q 

+0a«(z 2 P 2 - ziPi - Z)/3] , (A4) 

yafec ^gjjjg an antisymmetric tensor. The above color structure can be simplified by employing 
the identities, 

tr (T a T b T c ) = ^(d abc + if abc ) , d abc f abc = , f abc f abc = 24 , (A5) 

d abc being a symmetric tensor. 

In the collinear region with Z parallel to Pi, the terms proportional to g a p and g\ a do not 
contribute. Since the gamma matrices must be 7° = 7 + , 7^ = 7", and 7^ = 7_ = 7 + , the 
quark propagator between 7^ and 7^ is proportional to Zt, which is subleading in the collinear 
region. The factor g\ a indicates that 7 A must be 7", because of 7" = 7 + . According to 
Eq. (jHJ), the contribution from cZ(xiPi)7 A = d{x\P\)^{~ vanishes. The second term associated 
with gp\ in Eq. (|A4|) contains a collinear divergence. Due to 7" = 7 + , only the term — 2x 2 P 2a 
contributes. In the collinear region we have the approximation, 



2xoP 2n n_ 



1 



{x 1 P 1 - x 2 P 2 ) 2 (x 1 P 1 - a; 2 P 2 - If " n_ • Z LO1P1 - x 2 P 2 f (x 1 P 1 - x 2 P 2 - I) 2 J ' (A6) 
The expression in Eq. (IA3j) is then split, after the insertion of Eq. (jHJ), into two terms as 
shown in Eq. (|17|). 

The loop integral associated with Fig. 2(e) is given by 

x (f^^'^» ]^W • (A7) 

In the collinear region 7^, 7^, 7 and 7 a must be 7 + , 7", 7 T and 7 T , respectively. Using 
the eikonal approximation, 

^^JiMw-^^-r (A8) 
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and inserting Eq. @, Eq. (JSZI) leads to Eq. (ITHJ) . 

Following the similar treatment, the loop integral associated with Fig. 2(f) reduces to 

'l^^ iV^if • (A9) 

which then gives Eq. (fTDJ) . 

The integral of Fig. 2(g) is written as, 



= j ^d{x x P x ) la d{x 2 P 2 )u{x 2 P 2 )^ 



{p 2 -x 1 p 1 )^ p (p 2 -x 1 p 1 + iy^ 1 

(A10) 

In the collinear region with Z parallel to Pi, we have the sequence of the gamma matrices 
ryP _ ^ _ anc | ^ _ rj^g q Uar ]j propagator between 7 M and 7/3 gives the 
subleading contribution proportional to It- Hence, Fig. 2(g) does not contribute a collinear 
divergence. 

Figure 2(h) gives the loop integral, 



I [h) = -£h 7^A^h X 7^T^^d(x 2 P 2 )u(x 2 P 2 h 



2N C J (2vr)4 v A " (xxP-/) 2 ' v z/ v z/ ' (P 2 - x x P x y 
ti(T c T b T a )T c ^ a 

xyfefl) pw^lw^ • (A11) 

with the triple-gluon vertex, 

r rf« = /*»^ A (2Z - xxPi + x 2 P 2 ) a + g af} (2x x P x - 2x 2 P 2 - l) x 

+9x a (x 2 P 2 - X X P X - l) ] . (A12) 

In the collinear region with I parallel to P x , the terms proportional to g@\ and g\ a do not 
contribute. Because j x must be 7 + , the factor g@\ implies that 7^ must be 7" and the 
quark propagator between 7 A and 7^ is proportional to It, which is subleading. The factor 
gxa implies 7" = 7". Since 7^ is adjacent to the spinor u(x x P x ), it must be 7 + . The quark 
propagator between 7^ and 7" gives the subleading contribution proportional to It- Only 
the second term in Eq. (|A12|) contains a collinear divergence. 

Applying the approximation in Eq. (|A6|) . and inserting Eq. Q, Eq. (|A11|) becomes 

jw « E / ^itfi^i, ei)^^,^) , (Ais) 

O T J 



n=S,T 
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with the 0(a s ) collinear divergent PS piece, 



2m C F J (2vr) 4 v 1 1,1 (x 1 P 1 - I) 2 ' v 1 17 / 2 n_ • / 
<K£i - a:i) - 5 ^i - xi + p+ 



x 



(A14) 



The loop integral associated with Fig. 2(i) is written as 

T (i) ~eg 4 C F f dH - x x f x - / p x 2 P 2 - / j, , 

x " ( " 2P2) ^(f^fi^ v " (SlPl) !W^W ' (A15) 

After applying the eikonal approximation similar to Eq. (jA8|) . and inserting Eq. (|UJ), the 
above expression is simplified into 

j« « E / ^^eo^f (6,^) , (Ai6) 

n=S,T ^ 

with the collinear divergent PS piece, 

^ Wl) = 8^/ K (lll)7 (^-0 2 7 u^iPOp^^-orx). (A17) 
Figure 2(j) gives the loop integral, 



4iV c 7 (2tt)4 v 1 ^' (xiPx-/) 2 ' ^ " v * ^' a (z 2 P 2 -0 2 

x ^ (S-^o a 7 ^ iPi) p(^i - o 2 • (A18) 

Following the same procedure as for Eq. ()A15|) . we obtain 

« E / ^S(^,ei)^ 0) (ei^2) , (Ai9) 

n=S,T •* 

with the collinear divergent PS piece, 

At last, the integral associated with Fig. 2(k) is given by 

x (f^^ W.-U^ (A21) 

It is easy to find that a collinear divergence does not exist for the same reason as for Fig. 2(g): 
in the collinear region we have the gamma matrices 7" = 7+ and 7^ = 7 + , which is adjacent 
to the the spinor d{x 2 P2)- The contribution is then subleading because of the equations of 
motion for the final-state quarks. 
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APPENDIX B: 0(a s ) CORRECTIONS FROM FIG. 4 



In this Appendix we discuss the factorization of the initial-state collinear divergences in 
the 0(a s ) radiative corrections to Fig. 1(b), which are shown in Fig. 4. Figure 1(b) gives 
only the lowest-order PS distribution amplitude <f>g and the hard amplitude 

</4 0) (xi,£i) = — !— d(a;iPi)75u(xiPi)<J(^i - si) , 

4mn 



rr(o)/, „. x _ » 2 r m tr[ 7 ^(a: 2 P 2 )M(x 2 P 2 )7 /t (f 1 -x 2 AhV] mn 



because of 7 !/ (/fo+ /n~ — 1)7^ = 0, where the gamma matrices 7^ and 7^ come from the 
gluon vertices in Fig. 1(b). 

The integral from Fig. 4(b) is written as 

I [b) = ~\eg A Cl J -^dtxtPfoft *~ l )2 jJ(x 2 P 2 )u(x 2 P 2 ) 

Pi — X 2 P 2 Xi Pi+ / a ._ . 1 . . 

X7 ' (fi-W 7 (iTATTP 7 M(llA) /WW7ATv ' (B2) 

Following the same procedure as for Eq. (|A1|) . we obtain 

JW w Jdt^izu&Hjpfaxi), (B3) 



where the expression of 0^ (sci, £1) has been given in Eq. (|14|) . The factorization of Figs. 4(a) 

d 1 ) 



and 4(c) is performed in a similar way, leading to <p Sa Sc (xi, £1) in Eqs. (|T3Jl and (fT5|). 
The loop integral associated with Fig. 4(d) is given by 

rfdf) ^fl^ [ dH - x - nx P1-X2P2 a Pl+ / 

x^nix^)—— l p ; , (B4) 

r{xiPx - X2P2 - l) 2 {x 1 Pi - X2P2) 

with the triple-gluon vertex in Eq. ()A4|) . The same procedure as for Eq. (jA3|) leads to 



J« « J d^ d (x 1 ,^ 1 )Hl > \^,x 2 ) , (B5) 

with the function ^J^i, £1) shown in Eq. (fT7|) . 

The loop integral associated with Fig. 4(e) is written as 

x J^^ Va-W • (B6) 
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which is simplified into 

On the other hand, the loop integral associated with Fig. 4(g) is written as 



x ■ 



' M (Pi-x 2 P 2 ) 27/3 (P 1 -a; 2 P 2 + /) 27a 

xx ft + 1 



which reduces to 

,2 /• ^4 



1 8m N c J {2^Y d{XlFl)l \x l P 1 + lf 1 U[XlPl) P n„ ■ I (Px - x 2 P 2 + I)* " s 

This result differs from that of Fig. 2(g), which does not contain a collinear divergence. 

Neither Eq. (|B7|) nor (|B9J) is in the desired form. However, their combination is, as shown 

below, 

I (e) + I {9) ^ J^lfatJHMfaxi), (BIO) 

with the function <f>gl(xi,£i) shown in Eq. ((Hj) . 
The loop integral associated with Fig. 4(f), 



(/) eg'C F f dH - a x 2 P 2 + / - A-Xg P 2 

1 = ^VT 7 (27F d(XlPl)7 (x^ + O 2 7 ^^ 2 ^^^ 7 ^^-^) 2 7 " 



x i^^ tt( ^ ) P (x 1 A-U-/r (B11) 



leads to 



« | da0g(^i,6)4 o) (a^ 2 ) , (B12) 

where the function has been given in Eq. f)19|) . 

The integral associated with Fig. 4(h) is given by 

-teg 4 f dH -j, x x x P x -} j, * i — n \ Pi~x 2 P 2 



iih) = it I 



,m (pi-x 2 p 2 ) 2 

x 7 ^(x 1 P 1 ) ;2 , - ^— n2 7 Ap " — , (B13) 



tr(T c T 6 T a )r^ a Q 
'/ 2 (xiPi -x 2 P 2 -/) 2 (x 1 Pi-x 2 P 2 ) 2 ' 



with the triple-gluon vertex in Eq. (|A12|) . Following the same procedure as for Fig. 4(d), 
the integral becomes 

« | ^S(*i,£i)4 0) (^ 2 ) • (B14) 
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Similarly, we obtain 

J» « Jd^^&HP^xi), (B15) 

JW + /W« j d^]{x x , ^ x )Hf{^x 2 ). (B16) 

The functions and have been presented in Eqs. (jA14|) . 
(lATTf) and (jXlQj) respectively. 

Combining all the above contributions from Figs. 4(a)-4(k), we derive the 0(a s ) factor- 
ization in Eq. ()28|) . 



APPENDIX C: 0(a s ) CORRECTIONS TO B — > vr^ 

In this Appendix we discuss the factorization of the collinear divergences in Fig.4 for 
the semileptonic decay B — > 7rZP with the initial state and the final state being flipped. 
Figures 4(a)-4(c) can be factorized straightforwardly, leading to Eq. (|58|). Note that the 
hard amplitude from Fig. 4(b) depends on the loop momentum 

The loop integral from Fig. 4(d) is given by 

XI, tr(T<TW$ a 
*1 ^ Pflfc,-*^ + ')'(*. -W • < C1 > 
with the triple-gluon vertex, 

r% a a = f cba [g* P (2l + h- x 2 P 2 ) x + gp\(ki - x 2 P 2 - i) a 

+g Xa (2x 2 P 2 - 2ki - . (C2) 

The color factor is simplified according to Eq. (|A5|) . 

In the collinear region with I parallel to P 2 , only the term proportional to g a \ contributes 
by employing the argument the same as for Fig. 2(d) in Appendix A. Since 7^ must be 7"", 
only the plus component of k± survives. Applying the approximation similar to Eq. (|A6|) . 

-2k x p n +f3 



(h - x 2 P 2 ) 2 (h - x 2 P 2 + I) 2 n + ■ I 
and inserting Eq. ©, we obtain Eq. (fB5j). 



(h - x 2 p 2 y {h - x 2 p 2 + if 



(C3) 
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The collinear factorization of Fig. 4(e) can be achieved by applying the eikonal approxi- 
mation to the b quark propagator, 

K-M+I + m b _ 2(P 1 -fc 1 ) /3 - 7/3 (A- h-m h ) un 

npK p i-ki) « z-^—. b(Px-ki) 



(Pi -h + l) 2 - m\ IP y 1 17 2{P 1 -h)-l 

« -^r&(Pl - fcl) • (C4) 

n + • I 

The neglect of / is due to 7j g = 7 + in the collinear region. The second term on the right-hand 
side of the first line vanishes because of Eq. (|54j) . To derive the final expression, we have 
further dropped the power-suppressed terms proportional to k±. The integral associated 
with Fig. 4(e) then becomes 

He) f d A l 1- D \ a ^2 P2+ / 7/ D J ™+/3 



m N c J (2vr)4 v (x 2 P 2 + If '° y 2 2, Pn + -l 

v (P 2 - fci) 2 f-g* \ tr[ 7 "(P 2 - Ah^i^i ~ fci)rf(fc 1 ) 7a7 5 ] 

Note that the factor (P 2 — fci) 2 / (P2 — fci + Z) 2 in the second line indicates that Eq. ()C5j) has 

not yet reached the expected factorization form. 

The integral for Fig. 4(g) reduces, in a similar way, to 

T (g) _ W 2 f dH p x 2 P 2 + / - 1 n +l3 

1 ~ a aT / TTT^I u ( x 2 p 2)T f - p , n2 l5d{x 2 P 2 )-- 

8m N c J (27r) 4 (x 2 P 2 + I) 2 I 2 n + ■ I 

x (p 2 - £ + 2 l^" CFm v (p 2 - Wi - w ' (C6) 



Combining Eqs. ()C5|) and (|C6j) , we arrive at the desired factorization form in Eq. (J62J) with 
the collinear piece shown in Eq. (}64"|) . It is observed that the hard part is independent of 
the loop momentum /. 

The integral associated with Fig. 4(f) is factorized into Eq. (162)) straightforwardly with 
the collinear piece shown in Eq. ([65)1 . In this case the eikonal propagator arises from the 
approximation, 

( 1)7/3 (fcTTT) 2 w 2^-z = • (C7) 

Again, the neglect of / is due to 7j g = 7 + in the collinear region, and Eq. (}54|) has been 
employed to derive the final expression. 

The loop integral from Fig. 4(h) is written as 

<f S <W ,_ „. „ A- ft ,, n , >,„ , » 12 A- / 



(x 2 P 2 - Z) 1 



tr(T c T h T a )Tf a 



^-^'iw^ ' (C8) 
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with the triple-gluon vertex, 



r c i a Q = f cba [gpx(2l + h- x 2 P 2 ) a + g a p(ki - x 2 P 2 - l) x 

+gxa(2x 2 P 2 -2k 1 -l)p] . (C9) 

Following the same procedure as for Fig. 4(d), Eq. ()C8|) reduces to 

JW « | di 2 Hf{x u 6)0S(^2, 6) , (CIO) 

with the 0(a g ) collinear divergent function, 



/3 



X 



5(6 - £2) - 5 ( 6 - £2 + 



(Cll) 



The loop integral associated with Fig. 4(i) is given by 

"sl^^ptraF • (C12) 



Using the approximation similar to Eq. (|C7|) . the above expression is factorized as 

J« « | d^Vi, 6)4? (^2,62) , (C13) 
with the collinear divergent function, 

= 8^ / (S*<- P ^(tl^^4^ - ■ < C14) 

Following the same procedure as for Eqs. (|C5j) and (|C6j) . Figs. 4(j) and 4(k) lead to 
f dH „ (0) , . w _ D , x 2 /W „j, D \ 1 n +v (P 2 -A;i) 2 



Z 2 n+ • Z (P 2 - Jfei - J) 



Combining the expressions in Eq. (|C15J) . we arrive at the desired factorization form, 

I^ + I {k) ^ JdC 2 Hf\x 1 ,C2)4 ) j (x 2 ,^), (C16) 
with the collinear divergent piece, 
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We then discuss the factorization of the collinear divergences from Fig. 2 with the initial 
and final states being flipped, which represents the 0(a s ) corrections to Fig. 1(b). To 
simplify the discussion, we show only the PS parts below. The results are 



j(d) 



dH ( -g 2 C F \ trfUli ~ x 2 P 2 + /+ m fc ) 7 ^(P 1 - fc 1 )rf(fc 1 ) 7l /y 5 ] 

TTIq 



(2irf\ 2 "J [{Pi-x 2 P 2 + l) 2 -m 2 ](k 1 -x 2 P 2 f 



X 7T U ( X 2 P 2h {- p T m l5d(x 2 P 2 ) 

2m C F \x 2 "2 + I) l z n + ■ I 

2 r AA] x 2 P 2 + / 



2m C F J (2vr) 4 v z/ ' (:r 2 P 2 + If 

(C18) 



(e) /■ dH -g 2 C F \ trfn^n - x 2 P 2 + j[+m b )Yb(Pi - fciMfcih./? 5 ] 
™ J (2k)*\ 2 m °J [(p 1 _ X2 p 2 + /)2_ m 2 ](A;i _ ;r2 p 2)2 

X 8^N c U{X2P2h \ X2 P 2 + l) 2 l ^ X2P2) P^l ' (C19) 

« | d64 0) (^,6)0^(^2, 6) • (C20) 

Combining the above expressions, we have 

J (27r) 4 ^ 2 7 [(Pi-x 2 ^ + Z) 2 -m2](fci-a; 2 P 2 ) 2 

-ig 2 C F x 2 P 2 + / j. 1 n + 



4m (x 2 P 2 + l) 2 l 2 n + -l 



ig 2 C F f dH .., = n ^x 2 P 2 +/ 



u(x 2 p 2 h 



4m J (2tt) 4 v l > (x 2 P 2 + 1) 



Consider the loop integral associated with Fig. 2(g) under the approximation, 
P 1 -x 2 P 2 +/+m b Pi-x 2 P 2 + m b 2P 1/3 Pi + m b 



(C21) 



(Pi ~ x 2 P 2 + If - mf P {P 1 - x 2 P 2 f - m 2 ' (P, - x 2 P 2 + If - m 2 (P 1 - x 2 P 2 f - m 2 

2Pi • I n + Pi + m b 



(Pi ~ x 2 P 2 + I) 2 -m 2 n + -l(P 1 - x 2 P 2 ) 2 - mi 



The neglect of £ 2 /P 2 and x 2 /P 2 in the first and second propagators, respectively, is due to 
7^ = 7 + in the collinear region. The integral P 9 ^ reduces to 



/(<?) 



dH 2P 1 ■ I f-g 2 C F \ tr[^( A + m b )YHPi - h)d(h) lul 5 



(2tt) 4 (P x - :r 2 P 2 + If - m 2 \ 2 V J [(P ± - x 2 P 2 f - m 2 }(h - x 2 P 2 f 
-ig 2 C F p x 2 P 2 + / - 1 n + 

X ^m^ u{X2P2h (x 2 P 2 + lf l5d{X2P2) p—l ■ (C22) 
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The delicate combination of the first term in Eq. (jC21|) and Eq. (jC22|) leads to the correct 
factorization f 
to 5(£ 2 - x 2 ): 



factorization form with the hard amplitude Hg\x\, £2) and the collinear piece proportional 



f AC tr(0)/ c \ ~ l 9 C F f <Tl p X 2 P 2 + / 



L Tl-i. ' L 



4m J (2tt) 4 v ' *' 1 (x 2 P 2 + I) 2 

(C23) 



The sum of Eqs. (|H2T1) and (IH221) then gives Eq. (f7E|l . 
The loop integral from Fig. 2(h) is simplified to 



(2tt)^ 2 "J l( Pl - X2 p 2 )2_ m 2 ]{kl _ X2 p 2 + l y 



The loop integral associated with Fig. 2(i) reduces to 

J« « Jd^P(x u &)<l$(x 2 ,b), (C25) 
with the collinear divergent piece 4^sl{ x ^i £2) shown in Eq. (|C14J) . Figure 2(j) leads to 
jO) ~ / J!L ( Z9l CFm \ tr M p i -X2P2 + m b )j a b(P 1 - fei)d(A;i)7a7 51 



(2tt) 4 ^ 2 u y [(p 1 _ a ; 2 p 2 )2_ m 2 ](A;i _ a;2 p 2 + z) 2 

X 8^ M ^^^*^?^7 • (C26) 

Combining Eqs. (IU241) . (IH251) and (IH26L we arrive at 

„ j ^ (ll)6) {_g_ / |L u( , 2 p 2)75 |^ 7 ^ 2 p 2) ii^., (6 _ l2) 



(2tt)*\ 2 V J [(p 1 _ X2 p 2 )2_ m 2]( A;i _ X2 p 2 + / )2 

The integral associated with Fig. 2(k) becomes 
/(*) « / d4f 2Pl-/ ( zi- Czm^ ^[7^(A + ^)7 a W-^i)rf(^i)7«7 51 



(2tt) 4 (Pi - z 2 P 2 ) 2 - V 2 7 [(A - X2P2 + 2 - ™&] (h ~ x 2 P 2 + I) 2 

K^u(x 2 P 2 ) l5 ^ f 2 " • (C28) 



X 4m C F zy ,0 (x 2 P 2 - I) 2 ' ~ v ~" ZJ l 2 n + -l 
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Combining the second term in Eq. (jC27|) and Eq. (|C28|) . we obtain the correct factorization 
form in terms of the hard amplitude Hf\xi, £2) and the collinear piece proportional to 



The sum of Eqs. (1H271) and (IC281) then gives Eq. flZgJ). 
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Figure Captions 



1. Fig. 1: Lowest-order diagrams for tc-/* — > tc (B — > nlis), where the symbol x represents 
the virtual photon (weak decay) vertex. 

2. Fig. 2: 0(a s ) radiative corrections to Fig. 1(a). 

3. Fig. 3: (a)-(d) Infrared divergent diagrams factored out of Fig. 2(d)-2(k). (e) The 
graphic definition of the two-parton twist-3 pion distribution amplitudes. 

4. Fig. 4: 0(a s ) radiative corrections to Fig. 1(b). 

5. Fig. 5: (a) The Ward identity, (b) Factorization of 0(a^ +l ) diagrams as a result of 
(a). 

6. Fig. 6: Factorization of 0(a^ +l ) diagrams corresponding to Eq. (|3T)|) . 

7. Fig. 7: (a) A typical diagram of Gl^ 1 ^. (b) This diagram does not belong to G^n +1 \ 



43 



XiPi 



u 



X2P2 



d 



X1P1 



X2P2 



(a) 



FIG. 1 



1 




FIG. 2 



(b) (c) 



oo 

OO 



(e) 



FIG. 3 



3 




FIG. 4 



4 



"mum 
i 



+ 




+ 




(a) 



a 



N 



"my 



+ 



m5) 



a 



N 



(b) 



FIG. 5 



5 




6 




FIG. 7 



7 



